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Abstract 

The / — > d transition model by Duan and co-workers [Phys. Rev. B 66, 155108 (2002); J. Solid 
State Chem. 171, 299 (2003)] has been very useful in interpreting the / — ► d absorption, emission 
and nonradiative relaxation of light lanthanide ions in crystals. However, based on the assumption 
that the f N ~ l core spin-orbit interaction is weak compared to / — ► d exchange interaction, this 
model, in the original form, is not applicable to interpretation of the / — > d transitions of heavy 
lanthanide ions or actinide ions in crystals. In this work the model is extended to cover the cases 
of heavy lanthanides and actinides, where the spin-orbit interaction of / orbitals may be stronger 
than the f — d exchange interaction. 
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I. INTRODUCTION 



New luminescent phosphors for vacuum ultraviolet (VUV) excitation are required for 
plasma display panels and mercury-free fluorescent tubes, where the VUV emission from 
a noble gas xenon discharge is used to generate visible luminescence. Other applications 
where the VUV spectroscopy of lanthanides is involved are scintillator materials and VUV 
lasers. Due to these potential applications and availability of VUV excitation by synchrotron 
radiation, the VUV spectroscopy of lanthanide ions and actinide ions in crystal have recently 
become an important field of research. 

The VUV spectroscopy of lanthanide and actinide ions mainly involves the parity allowed 
nf N <-» n/ iV_1 (n + l)d transitions (n = 4 for lanthanides and 5 for actinides), including 
ground and excited state excitations and emissions. It is well known that the 4/ — > 4f N 
spectra of lanthanide ions in noncentrosymmetric environments are dominated by many 
sharp zero-phonon lines. Interpretation of the energy levels and transition intensities may 

n nn 

be modeled with a empirical crystal-field Hamiltonian [1| and Judd-Ofelt theory|2j,y] respec- 
tively. By contrast, the spectra of transitions between the 4/^ and 4f N ~ l 5d configurations 
are feathered by broad-band structures, with some zero-phonon lines only resolvable at 
temperature below liquid nitrogen. This is due to great difference in vibrational equilib- 
rium between 4f N ~ l 5d and 4/^ configurations, which makes the the transition intensities 
be dominant by vibronic broad bands. The phenomenological crystal-field Hamiltonian for 
Af N configuration has been extended to 4f N ~ 1 5d configuration by adding the crystal-field 
and spin-orbit interactions for 5d orbitals, and Coulomb interactions between 4/ Ar_1 core and 
5d orbitalsj^]. The 4f N <-> 4f N ~ 1 5d transitions are electric dipole allowed, whose relative 
rates can be calculated straightforward. Actually extensive calculations have been carried 
out for trivalent lanthanide ions in Cry8 ta ls flfi which give 8 at lsfaC to ry a„t with 
experimental 4/ — > 5d spectra. However, the calculations give hundreds to thousands tran- 
sition lines which convolve into several to a few dozen of broad bands after taking vibronic 
transitions into account to wash out the fine structures, making straightforward interpreta- 
tions of the simulations and straightforward predictions of changing of broad-band numbers, 
their positions and relative intensities with ions and crystals quite difficult. 

Recently Duan and co-workers 0,0,0] simplified the calculations by considering only the 
main interactions in the 4f N ~ l 5d configuration. The model gave quantum numbers char- 
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acterizing various groups of states, transition selection rules, and an expression of relative 
intensities with these quantum numbers. Application of the model to Eu 2+ and Sm 3+ in 
various crystals successfully explained the measured spectra. However, in the model there 
is an implicit assumption that the exchange Coulomb interaction between / and d orbitals 
is stronger than the spin-orbit interaction in the f N ^ core. This assumption no long holds 
for heavy lanthanide ions or actinide ions in crystals. 

In the present work, the original model is extended to the case where f — d exchange 
Coulomb interaction may be weaker than that of the f N ~ x core spin-orbit interaction. The 
energy and transition line strengths for / — d transitions directly applicable to actinides and 
heavy lanthanides are given analytically. 

II. EIGENFUNCTIONS AND EIGENVALUES 

A. f N configuration 

The f N crystal-field splitting is well-understood both experimentally and theoretically 
via phenomenological crystal-field simulation. However, in f N f N ~ 1 d transitions, the f N 
crystal-field energy levels are either hard to distinguish from vibronic lines in low temperature 
or unresolvable at all at temperature higher than 100K due to vibronic bands. In this work 
we are aimed to give an model which may interpret number of broad bands, their positions 
and relative intensities. We neglect f N Crystal-field splitting. In such a case, the highly 
degenerate eigenstates can be written as 

| [r}SL]J) = C \rjSLJ) + £ a I^L, J) . (1) 

i>l 

Usually the states are denoted with the label of the main component \i]SLJ), whose coeffi- 
cient is usually close to 1 for lanthanide ions. 

B. f N ~ 1 d configuration 

The general interactions for the f N ~ l d configuration can be written as 

H(f N ^d) = H Coul (ff) + H c{ (d) + H Coul (fd) + H so (f) + H so (d) + if d (/). (2) 
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The first two terms are the strongest terms of the following form 



Ho(f N ~ 1 d)= £ F k (ff) £ C k (i).C k (j)+Y: E B k q Cf\d), (3) 

fe=2,4,6 l<i<j<N-l k=2,4-k<q<k 



where Fk(ff) are slater integrals, which are usually treated as adjustable parameters. B k 
are crystal-field parameters for the d electron. Note that only those B h with (k, q) allowed 
by the site symmetry are nonzero. Ho(f N ~ 1 d) contribute to splitting of f N ~ x core into 
energy levels characterized with spin and orbit angular momenta S and L, and d orbital into 
strong crystal-field energy levels characterized with site-symmetry labels. The contribution 
to splitting from the third term of (J2J), i.e., the Coulomb interaction between J^ -1 core and 
d orbitals, can be approximate with an isotropic exchange term 

H exc (f N ~ 1 d) = -J C xcSf • S d , where 

Jexc = ^GiUd) + ^G 3 (fd) + ^G 5 (fd). (4) 

Here £71,3,5 (/d) are / — d Coulomb exchange radial integrals. It is straightforward to show 
that if exc commutes with the total spin of f d. 

The following approximation is often used for /-electron spin-orbit interaction within a 
given zero-order J^ -1 core energy level characterized by rjSL: 

JV-l 

H so (f) = E ^ nlSi ' 11 ~ KslSi ■ Li, where (5) 
i=i 



A 



rjSL 



/(/ + 1)(2/ + 1) , 7 iV„.aniT,Hii,iV 



\ S(S + l)(2S + l)L(L + 1)(2L + r 



{raSLWV^WraSL)^- (6) 



In the case that Sf takes the largest possible value for the f N ~ x configuration, X Vf s f L f 
is simply sign(8 — N)£f/2Sf Q, IJ- In general, H so (f) commutes with the total angular 
momentum of J^ -1 core, no matter the above approximation in Eq. (J3J) is used or not. 

Other terms are not important in the interpretation of broad bands in / — d spectra and 
neglected. Therefore we have the effective Hamiltonian written as the sum of the above 
important terms as 

H cS = Hoif^d) - H exc (f N ' l d) + #«(/)■ (7) 

It is straightforward to check that H e g commutes with the following effective "angular mo- 
mentum" operator 

J cfr = S / + S (i + L / . (8) 



Note that this operator is not the total angular momentum operator for f n ~ l d, since it does 
not contain the orbit angular momentum of d electron which is usually quenched in low 
symmetry sites. 

Former work used an implicit assumption that fd exchange interaction is stronger than 
f N ~ l spin-orbit interaction. In that case the energy levels were written as 

E (jl^rjSfLf, 2 d>; Sj)) = Eo(vS f L f ) + t di - J cxc \s{S + 1 - S f (S f + 1) - ^ 

( 2S + 1 \ J{J+l)-S(S + l)-L f (L f + l) 
+ { 2 ~ 257TT ) **' 8 ' L * 2 ' (9) 

However, the strength of the exchange interaction decreases as the nucleus charge increase, 
while at the same time the spin-orbit interaction increases. It happens that for heavy 
lanthanide and actinide ions, the cases where spin-orbit interaction is comparable or even 
stronger than exchange interaction need to be considered. In the case H so (f) is much 
stronger than H exc (fd), opposite to the one considered by Duan et. aZp]], The coupling 
(SO^ f x ) x S*Of d may be preferred and the approximate eigenstates can be written as 
(l N ~ l rjSfLf, Jf), 2 df, J^> and the eigenvalues can be written as 

E( v S,L f J„ %; J) = E ( V S,L„ %) + \, SfL { J ' { J < + S > {S > + L ' {L > + l) 



_ ^ JM + i) + g^ ++ i)-ML / + i) [J(J + 1) _ JAh + 1} _ Sd(Sd + 1}] (10) 

In the medium case where H so (f) and H exc (d) are comparable, The effective of H e ^ can be 

l N ^ 1 rjS fLfJf, 2 di J) or ^^rjSfLf, 2 diSj\&s bases. Eigenvalue and 



calculated with either 



wave-function of each eigenstate can then be obtained by diagonalizing the matrix. Here we 
give the matrix element of effective Hamiltonian under the bases l N ~ 1 r]SfLf 2 diSJ) as 



r- l r)S f L f , %;SJ H cS r- l r]S f L f , 'd^S'J 

x rr x £ \Cc ItN-l a r \ , \ t S(S + 1) - Sf(Sf + 1) - S d (s d + 1) 

= S diA {5s fS ' f h f L' f Sss'[(Eo(f vS f L f )+e dt )-J cxc ^ '- ^ '- * '-] 

+ U~n J+Lf+S ' f )} + + ±)(vS f L f SJ\\V u \\ V 'S' f L' f S'J) (11) 

Under the approximation in Eq.fjSJ), the matrix for H e R reduces into many 2x2 blocks and 
the diagonization become straightforward. 
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III. ONE-PHOTON TRANSITION LINE STRENGTH BETWEEN f N AND f N - l d 

The / to f N ~ 1 d transitions are electric dipole allowed. Here we consider only this 
mechanism. The electric dipole moment is a spin independent rank 1 tensor in both total 
orbital angular momentum and total angular momentum spaces. It can be written as 

£ = £^f> ? (i), (12) 

9 i=l 

where e q is the q component of the polarization vector and r q (i) is the q component of the 
position of i th electron. 

Using the second quantization techniques, the electric dipole momentum can be written 

as 

D (o.i)i = y /2(f\r \d) {[(a+fr^d^}^ - [( a +)(V2-2) 5 (V2-3)](o-i)i 3 } , (13) 

where 

and (a+)( srasim ') are components of tensors that transform under symmetry operator the 
same way as basis \smjmi), and (f\r\d) is radial integral. The coupling of two creation- 
annihilation operators are just coupling of two tensors to give a new tensor. In Eq.(|13|). a + 
and a couples to give a rank (0 • 1) tensor of spin and orbital angular moments. 

Using coupling and recoupling techniques, we can rewrite the electric dipole momentum 
into the following two forms. 

D = E Ci g2 [(a + )( 1 / 2 - 3 )s( 1 / 2 - 2 )](°- 3 - 2 )^' 2 ^) + ■ ■ • (15) 
= E C 2 fgi9 J(a + )( 1 / 2 - 3 )a( 1 / 2 - 2 )]^- 1 / 2 - 2 )^ 1 ' 292 ) + ---, (16) 

Of, 91 ,92 

where the neglected terms (• • •) will not contribute when / states are on the left and f N ~ 1 d 
are on the right and hence are not written ou t explicitly. C qiq2 and C 2 /gi(J2 are appropriate 
coefficients that depend on (qx, q 2 ) and (jf, qx, q 2 ) respectively. The matrix elements of 



D between initial states 



f N d°SLj) and final states {f N - 1 d 1 {{S f s d )SL f )J' 1 d l ) ( or final 



states 



(/ 1 d l ((SfLj)JfS ( i)J\dij) can then be obtained via Wigner-Ekwart theorem 



0. 
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Transition line strength S between initial states \Ii) and \Ff), where i and / are indexes 
to distinguish degenerate states, are defined as follows 



s(i,F) = j:mD\Ff}\ 2 . 

i,f 



(17) 



Using orthonormal relations of coupling and recoupling coefficients [n|, after a lengthy but 
straightforward analytic calculation, we finally get the line strength for isotropic absorption 
or emission. For the case I of strong H exc (f N ~ 1 d) which is applicable to light lanthanides, 
the line strength is: 



S iso (\f N V SLj) <- [[{f^TjfSfLf, 2 d i ;S'J' 



35 I j> j s' 



(18) 



where the [di\ is the degeneracy of di crystal field levels, and [S] etc. are short for (25* + 1) 
etc., respectively. 

For the case II of stronger H so (f) than H eyiC (f N ~ 1 d), the line strength is: 



f N V SL.j)^ [{f N - x ri,SfLf)J f , 2 d l ;J f ))=N(f\r\d) 2 [d i ][S,L,J f ,J,.r 



x(r v sL{\r-\s f L f ) 



t 



7/2 

E (-i)^ 7/2 b/] 

jf =5/2 



V 




\ 



,(19) 



IV. CONCLUSION 

The model for f N — > f N ~ 1 d transitions proposed earlier|7j has been extended to the case 
II where spin-orbit interaction in f N_1 is stronger than the isotropic exchange interaction 
between f N_1 and d by utilizing Racah-Wigner algebra and second quantization techniques. 
The result is expected to be useful for actinide ions where the effect of / spin-orbit interaction 
is stronger than the exchange interaction. Heavy lanthanides fall into the medium case where 
the effect of spin-orbit interaction is compariable or slightly stronger than the exchange 
interaction, where the case II may serve as an approximaition. 
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